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o : 

, We study transport through an electronic Mach-Zehnder interferometer recently devised at the 

, ^ ■ Weizmann Institute. We show that this device can be used to probe statistics of quasiparticles in 

the fractional quantum Hall regime. We calculate the tunneling current through the interferometer 
, as the function of the Aharonov-Bohm flux, temperature and voltage bias, and demonstrate that its 

flux-dependent component is strongly sensitive to the statistics of tunneling quasiparticles. More 
' specifically, the flux-dependent and flux-independent contributions to the current are related by a 

t—\ , power law, the exponent being a function of the quasiparticle statistics. 

I. INTRODUCTION 

i One of the key features of the quantum Hall effect (QHE) is the fractional charge and statistics of quasiparticles. 
The seminal shot noise experiments of the Weizmann and Saclay groups 1 '^ allowed for direct observations of fractional 
charges. A recent experiment on the mutual fractional statistics of two quasiparticles with different charges has been 
published in Ref. 0- It involved a setup consisting of an island of a fractional quantum Hall liquid embedded in a 
liquid with a different filling factoid. Theoretically, several approaches'*^ have been proposed as possible tools to 
probe mutual statistics of identical quasiparticles. However, to this date, no experimental verification of the statistics 
I \ of identical quasiparticles has been reported. 

In this paper we propose a different approach to observing the statistics of identical quasiparticles. Our approach 
employs the electronic analogue of the Mach-Zehnder interferometer, recently designed at the Weizmann Instituted. 
This device has been used to observe the Aharonov-Bohm effect in the integer quantum Hall regime 8 . As is shown 
below, in higher magnetic fields this type of device would allow the observation of fractional statistics. 

The proposed method of the observation of fractional statistics has a number of advantages in comparison with 
other set-ups. In Refs. 00 one needs to measure the noise or the current-current correlation function. In our approach 
it is sufficient to find the current through the interferometer. This resembles Ref. 0. However, in Ref. one has to 
control the number of quasiparticles trapped in the interferometer in order to probe fractional statistics. There is no 
(f-) \ such difficulty in our case. Besides, the interference pattern can be observed in the Mach-Zehnder interferometer at 
\Q . a larger interferometer size than in the standard geometry'. 

An electronic Mach-Zehnder interferometer is sketched in Fig. 1. Charge propagates along two quantum Hall edges 
and tunnels between the edges at the point contacts QPC1 and QPC2. Figs, la) and lb) depict two possible setups: 
a) tunneling takes place between two different fractional quantum Hall puddles; and b) tunneling is between the edges 
of a single puddle. In the latter case, fractionally charged quasiparticles can tunnel at QPC1 and QPC2. In the 
former case, only electrons are allowed to tunnel 9 . Pinching off the edges closer to each other, one can deform the 
system depicted in Fig. lb) into the configuration Fig. la). The tunneling current between the edges depends on 
. their voltage difference and the magnetic flux through the region A-QPC1-B-QPC2-A. As shown below, the current 
includes a flux-independent contribution Iq and a contribution /$ which oscillates as function of the magnetic flux 
with period <&o = hc/e. We will calculate /$ and Jo as functions of tunneling amplitudes Ti and T2 at QPC1 and 
O QPC2. If no tunneling occurs at QPC2, i.e. T2 = 0, then the total current I(Ti, T2) = Io(Ti, 0) is independent of the 
\ magnetic flux. At weak tunneling at QPC2, i.e. small T2 < Ti, the flux-dependent component of the current scales 
as I$(Fi, T2) ~ [Io(ri,r2) — Iq(Ti, 0)] b , where the exponent b depends on the statistics of tunneling particles. In the 
case of electron tunneling, b = 1/2. In the case of the tunneling of fractional quasiparticles, the exponent depends on 
• their statistics and is always greater than for electrons. 

The outline of the paper is the following. First, we give a qualitative explanation of our results. Next, we introduce 
the model which is used in our calculations. Then, we calculate the I — V curves, respectively, for quasiparticle 
tunneling in the fractional QHE regime with filling factors v = \/{2m + 1), and electron tunneling between quantum 
Hall liquids with filling factors ^2 = l/(2mi 2 + 1). In the latter case our calculations follow the standard rout o 5 i 10 . 
In the former case one has to carefully treat the Klein factors describing quasiparticle statistics. The Appendices 
discuss some technical details of the perturbation theory employed. 
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II. QUALITATIVE DISCUSSION 

We study charge transport through the Mach-Zehnder interferometer in the limit of weak tunneling at QPC1 and 
QPC2. Thus, we assume that the edges of the fractional quantum Hall liquid are far from each other in comparison 
with the magnetic length l . Hence, r 1; T 2 <C ?i/t c , where the ultra-violet cut-off scale r c ~ l /v and v is the excitation 
velocity along the edges. The tunneling amplitudes can be controlled by gate voltages. We assume that the gate 
voltages are chosen such that one of the situations shown in Figs, la) and lb) takes place. Other situations, e.g. weak 
quasiparticle tunneling at QPC1 and weak electron tunneling at QPC2, are also possible but will not be considered 
here. 

We first address the electron tunneling case depicted in Fig. la). There are three contributions to the current 
I = 1% + 1% + Ii2- The first contribution is due to the tunneling at QPC1, I± = ci |Ti | 2 , where ci depends on the 
temperature and voltage and is independent of the magnetic flux through the interferometer. The second contribution 
is due to the second quantum point contact, I 2 = ci |T2 1 2 - Finally, a contribution arises due to quantum interference 
between the electrons which tunnel at QPC1 and QPC2. This contribution equals I\2 = C2riT2 exp(— icf>) + c.c, where 
c 2 depends on the voltage, temperature and the distance between the point contacts, and </> is the difference of the 
Aharonov-Bohm phases picked up by electrons propagating between QPC1 and QPC2 along two edges. The phase 
= Ic[Jqpci-A-qpc2 — Jqpci-b-qpc2 = 27r<3> /^o, where A is the vector potential, $ < the magnetic 
flux through the region A-QPC2-B-QPC1-A, e < the electron charge, and $0 the magnetic flux quantum hc/e. 
The coefficients c\ and C2 will be calculated in section IV. Thus, the current exhibits a periodic dependence on the 
magnetic flux $ with period $0 ■ The amplitude of the flux-dependent contribution to the current = max J12 is 
related to the flux-independent contribution Iq = I\ + I2 by the equation 

/*(ri,r 2 )~ [io(ri,r 2 )-/ (ri,o)] 1 / 2 . (l) 

Eq. JU is derived rigorously in section IV. 

A similar relation with a different exponent can be obtained for the set-up Fig. lb) in which fractionally charged 
quasiparticles tunnel between the edges. Quasiparticles in a quantum Hall liquid with the filling factor v — l/(2m + l) 
can be described as point charges q = ve with attached solenoids^. Each solenoid carries one magnetic flux quantum 
$o- When one quasiparticle makes a circle around another it picks up the Aharonov-Bohm phase 

9 = 2ttv (2) 

which describes fractional statistics. The total flux <t through the interferometer includes the contribution from the 
applied magnetic field 3? = J BdS and the statistical contribution from the flux tubes attached to quasiparticles. 

Fig. 2 is topologically equivalent to Fig. lb). As is clear from Fig. 2, the solenoids attached to quasiparticles 
do not contribute to the magnetic flux through the dashed circle, if there is no tunneling between edges lq and 2q. 
Thus, in the absence of tunneling the total magnetic flux $ = $. Each tunneling event changes the flux $ by one 
flux quantum. The flux decreases by |$o| when a quasiparticle tunnels from the outer edge to the internal edge, and 
increases by |$o| for the tunneling events from edge 2q to edge lq. Hence, $ = $ + n<i>o, where n is an integer. 

We will assume that the electrochemical potential of source SI is higher than the electrochemical potential of source 
S2 by eV, where e is an electron charge. These electrochemical potentials are equal to the chemical potentials of edges 
lq and 2q. 

Let us first consider the simplest limit of zero temperature. In this case, quasiparticles can tunnel from the edge 
with the higher chemical potential to the edge with the lower chemical potential only and cannot tunnel from edge 
2q to edge lq. The tunneling probability can be derived in the same way as Eq. Q , 

p($) = SidT!! 2 + |r 2 | 2 ) + (c 2 r 1 r;ex P (-^) + c.c), (3) 

where <j> — 2iri/Q/$o = 2ttv&/&o + 2mm is the Aharonov-Bohm phase accumulated by a quasiparticle with the 
charge ve along the path A-QPC2-B-QPC1-A. The phase cf) is periodic in n with period 1/v. Let us assume that 
initially n = k/v. Then p($) = p(&) and the transfer of a quasiparticle requires the time ti = l/p(<i>). After a 
quasiparticle tunneling event the statistical flux changes and the second quasiparticle tunnels after the time interval 
i2 = l/p( ( E > + $o)- After 1/v tunneling events we return to the situation with n being a multiple of \ jv. The time 
needed for the transfer of \ jv quasiparticles, i.e. a single electron charge, is t ~ X^=i h- Hence, the current 
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where Ii — vep($ + [I — l]$o)- If the quasiparticles did not carry flux tubes the current through the Mach-Zehnder 
interferometer in the presence of the magnetic flux ($ + [I — l]$o) would equal Thus, the current of quasiparti- 
cles obeying fractional statistics is a harmonic average of l/v currents corresponding to the imaginary situation of 
fractionally charged quasiparticles which do not obey fractional statistics. 

If r 2 <C Ti then the current / can be expanded in powers of r 2 . The term of order T^T™ is proportional to 
exp(27wi[u — «;] $/$())■ The current Q is a periodic function of the magnetic flux with period $o- Hence, [u— w] = k/v 
in all non-zero contributions. Thus, the linear in | T2 1 contribution to the current vanishes. Hence, the leading power 
of r 2 in the expansion is | T2 1 2 - Such second order terms do not exhibit a magnetic flux dependence. Hence, the 
flux-independent contribution to the current Iq satisfies the equation 

/o(r l5 r 2 ) -/ (ri,o) ~ |r 2 | 2 . (5) 

The leading flux-dependent contribution scales as 

/ # ~ consti^/" exp(27ri$/$ ) + c.c. (6) 
Finally, one gets from the comparison of Eqs. (J5J and JBJ at r 2 <C T\ 



j*(Ti,r 2 ) ~ [/o(ri,r 2 )-/ (ri,or /e , (7) 

where 9 is the statistical angle {J2J. Eq. Q is derived rigorously in section V. The flux-dependent and flux-independent 
contributions to the current can be expressed via the maximal and minimal values of the current as function of the 
magnetic flux: la = [max/ + minJ]/2, /$ = [max/ — min/]/2. The exponent ir/9 is determined by the quasiparticle 
statistics. It equals 1/2 for fermions, Eq. JIJ, an d exceeds 1/2 for fractional quasiparticles. 

Luttinger liquid effects are know to give rise to power laws in the edge physics in quantum Hall systems^. We 
would like to emphasize that Luttinger liquid physics is irrelevant for the above result, Eq. Q. The exponent n/0 
emerges due to the statistical magnetic flux carried by quasiparticles, i.e. due to their fractional statistics. 

We now discuss the finite temperature regime. In this case quasiparticles can tunnel both from edge lq to edge 2q 
and from edge 2q to edge lq. The tunneling probabilities are related by the principle of detailed balance 



+ [k x - + r]$ ) = 7P+($ + [kx - + r- 1]$ ), (8) 
v v 

where 7 = exp(— i/eV/kgT); 1 < r < v; the convention r — 1 = 1/v is used at r — 1; p+(3> = $ + 
statistical flux before tunneling) and = $ + statistical flux before tunneling) denote the probabil- 

ities of the tunneling events from edge lq to 2q and from edge 2q to edge lq respectively. They depend quadratically 
on I\ 2 and are calculated in section V. The tunneling probabilities depend on the total magnetic flux which includes 
the statistical contribution $ s = [k x i + s]<i>o, 1 < s < l/v. Since the Aharonov-Bohm phase due to the statistical 
flux Q s equals 4>ab = (27r/c + 2-7rz/s), the probabilities depend on $ and s only and are independent of k. Thus, the 
tunneling current is given by the equation 



i/ v 

I = ev^fr ]p+(* + r$o) - P- ($ + r* )] , (9) 

r=l 

where f r is the probability to find the system in one of the states with $ s = [k x i + r]$0j 1 < r < > k being an 
arbitrary integer. The distribution function f r can be determined from the steady state condition 



/ r [p+(* + r$ ) + r$ )] - /r-U»+(* + [r - l]«o) + fr+iP-i® + [r + l]*o)], 



(10) 
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where we use the convention 1/v + 1 = 1. Using Eq. (JSJ) the kinetic equation (|10fl can be rewritten as 



+ r$„)[/r - 7/r+l] = P+(* + [r - l]*o)[/r-i - jfr]',r = 1, . . . , 1/V 
One finds from the above system of equations 



(11) 



fr - 7/r+l 



p+($ + r$ )' 



(12) 



where a is a constant independent of r. To calculate this constant we add Eqs. 1)12(1 with all possible r. Since 
Y,fr = 1, one gets 



Thus, the current, Eq. ©, equals 
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Z^r=l p + ($+r<I>o) 



(13) 



ej/ X! p +(* + r$ o)[/r - 7/r+i] = ea 
r=l 



1/1/ 



Z^r=l I' 



(14) 



where I' r = ue(l — 7)p+($ + r$o)- Similarly to Eq. (0J, the currents I' r can be understood as the currents of 
fictitious fractionally charged quasiparticles which do not obey fractional statistics. A rigorous derivation of Eq. I|14|) 
is discussed in section V. Finally, the same analysis as in the zero-temperature case shows that Eq. Q is satisfied at 
finite temperatures. 

Eq. is the main result of the article. An experimental test of this relation will allow the observation of fractional 
statistics. In such experiment one needs to change T2 at fixed T\. The tunneling amplitudes are controlled 
by gate voltages. Generally, any change of gate voltages affects both tunneling amplitudes. They can be controlled 
independently only if QPC1 and QPC2 are far from each other. In most interferometer set-ups an increase in the 
distance between QPC1 and QPC2 would result in the suppression of the interference pattern. Fortunately, this is 
not the case for the Mach-Zehnder interferometer. The calculations for the electronic Mach-Zchnder interferometer in 
the integer quantum Hall regim e 13 ! 14 show that the visibility of the interference pattern depends not on the distance 
between QPC1 and QPC2 but only on the difference of the distances between the point contacts along two edges. We 
confirm the same conclusion for the fractional quantum Hall systems in sections IV and V. 

In order to determine the relation between the flux-dependent and flux-independent components of the current 
one has to vary the magnetic field and measure the current at different values of the field. The magnetic flux $ 
through the region QPCTA-QPC2-B-QPC1, Fig. lb), includes two contribution: the flux through the hole in the 
interferometer, i.e. the upper half of the region QPC1-A-QPC2-B-QPC1, and the flux through the lower half of 
the region QPC1-A-QPC2-B-QPC1 which is occupied by a quantum Hall liquid. If the magnetic field changes at 
the fixed density then the filling factor in the lower half deviates from v. As a result, quasiparticles can enter the 
region QPC1-A-QPC2-B-QPC1. Each of them brings one flux quantum. This does not change any of the results 
of the paper. Indeed, we predict that the current is a periodic function of the magnetic flux through the region 
QPC1-A-QPC2-B-QPC1 with period $o- Hence, changing the flux by one flux quantum does not affect the current. 

Several other quantum Hall interferometer set-ups have been discussed in the literature. The simplest set-up^ 
is illustrated in Fig. 3. In contrast to the Mach-Zehnder interferometer, the effective magnetic flux perceived by 
quasiparticles does not change after tunneling events in the set-up Fig. 3. Hence, the current is independent of the 
statistical phase 8, Eq. J2J (if n0 quasiparticles are trapped between the quantum point contacts). On the other 
hand, the current exhibits a "fractional" Aharonov-Bohm periodicity with period $0/^- On the technical level the 
set-up^ and our problem are described by very similar models (see section III). The main difference consists in the 
Klein factors which describe fractional statistics. They are present in the model of the Mach-Zehnder interferometer 
(section III) and are absent in the model of Ref. Of This difference between the models results in qualitatively 
different transport behavior. A set-up related to the Mach-Zehnder interferometer was studied in Ref. 0, Fig. 4. 
The interferometer, Fig. 4, has the same topology as the Mach-Zehnder interferometer but includes three edges and 
three quantum point contacts. In the absence of QPC3 the set-up Fig. 4 is equivalent to Fig. 3. If weak tunneling 
at QPC3 is allowed the system exhibits strong telegraphic noise which carries information about fractional statistics. 
In contrast to Ref. we investigate the set-up with two point contacts that was studied experimentally^. In our case 
not only the noise but also the average current is strongly sensitive to fractional statistics. 
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III. MODEL 



A. Effective action for electron tunneling case 



We first consider the system represented by Fig. la). Its low-energy behavior can be described by the chiral 
Luttinger liquid modelS.. Fig. 5 illustrates the model. The two edges 1 and 2 correspond to two chiral Luttinger 
liquids with the same propagation direction. The dashed lines describe quantum point contacts where electrons tunnel 
between the edges. Note that the respective distances L and L + a between the point contacts along the two edges 
are different. The Lagrangian assumes the standard form 9 

L = f dxdt { d t<t>kd x <t> k + v{d x 4> k ) 2 ] - f dt{T x + T a ), (15) 

where T% and T 2 are tunneling operators, v is the excitation velocity along the edges, and {</>&} represent two chiral 
Bose fields which satisfy the following commutation relations 

[4>i{xi,t = 0), (f>p(x p , t = 0)] = iir5i p sigii(xi - x p ). (16) 
The Bose fields are related to the charge densities pi through 



Pi = Wuie/2ir)d x (tn, (17) 

where e is an electron charge and vi = l/(2m; + 1) are the filling factors of the QHE puddle defined by edge 2e and 
the QHE strip bounded by edge le. The tunneling operators^ 



Ti = riexp(t[0i(O,i)/VEl - M^,t)/V^\) + h.c; 
T 2 =T 2 eMi[ML,t)/^~ML + a,t)/^}) + h.c (18) 

are proportional to the electron annihilation and creation operators ~ exp(±i<^; / \Jv~x) ■ The action also includes 
operators describing simultaneous tunneling of several electrons but they do not play a significant role at low energies. 

In the presence of a magnetic flux $ through the region A-QPC2-B-QPC1-A, the tunneling amplitude T 2 should 
be multiplied by the phase factor exp(27ri$/$ )j where $0 = hc/e is the flux quantum™. This phase factor describes 
the difference of the phases e/(Kc) Jq PC1 Adl accumulated by the electrons moving along two edges between the 
point contacts. For one edge the integration path is QPC1-A-QPC2, and for the other edge the integration path is 
QPC1-B-QPC2 (Fig. 1). 

The voltage bias V results in the chemical potential difference fi\ — ^ = eV between edges le and 2e. We would 
like to emphasize that the potential difference between the edges is determined by the voltage difference between the 
sources SI and S2 in both set-ups Fig. la) and Fig. lb); however, the shape of the edges le and 2e, Fig. la), is 
different from the shape of the edges lq and 2q, Fig. lb). We will use the interaction representation which makes 
both chemical potentials equal and introduces time-dependence into the tunneling amplitudes: 



r 1 ,r 2 ~ex P ( — — ). (19) 

The difference of the chemical potentials results also in the difference of the average densities between the edges. We 
assume that edge 2 is connected to the ground. Then the average density p 2 ~ (d x <t>2) is zero in the limit of weak 
tunneling. It will be convienient for us to shift the Bose field on the first edge 4>i(x) — ► <fii(x) — f(x) by a function 
f(x) of the coordinate in such a way that (d x (f>x) becomes zero. The mean charge density on edge 1 is proportional 
to the voltage bias and can be found from the minimization of the Hamiltonian 



H = h dx 



in 



{d x ct>i? 



2wh 



(20) 



This yields q — (d x 4>i) — eV-JU\/(hv). We next shift the field <fii — > 4>i — qx such that (d x (j>i) vanishes. At the same 
time T 2 is multiplied by the factor exp(ieVL/hv). 

Hereafter we assume that v\ > v 2 . The case v\ = v 2 = 1 corresponds to the free electron probler 



,13.14 
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B. Effective action for quasiparticle tunneling case 



We now consider the system depicted in Fig. lb). The model is given by the chiral Luttinger liquid action H15|) 
with a different choice of the tunneling operators^. Several modifications immediately follow from the fact that the 
quasiparticle charge ve differs from the electron charge. 1) The tunneling operators should be expressed via the 
quasiparticle annihilation and creation operators ~ exp(±iy^(pi); 2) the flux-dependent phase factor in T 2 is now 
exp(27r^i$/<f>o); 3) the phase factor exp(ievVL/frv) should be used instead of exp(ieV L /hv); and 4) in the interaction 
representation the time-dependence of the tunneling amplitudes becomes 



, ievVt „ 

r 1 ,r 2 ~ex P ( — — ). (21) 

The fifth difference from Eqs. Ijl8(l consists in the introduction of two Klein factors k\ and k 2 (a related model without 
Klein factors has been considered in Ref. Ilfifl : 



T? = Tmi exp{iy/v[4n(0, t) - 0a(O,i)]) + h.c; 
T§ = T 2 K2 exp(iy/v[<tn(L, t) — fa(L + a, t)]) + h.c, (22) 

where the commutation relations are 



k±k2 = exp(— 2ttvi)k2Ki; 

K1K2 = k-2 nx exp(27rfi). (23) 

These commutation relations can be understood from a locality argument similar to Ref. 7. Two tunneling 
operators affect two distant parts of the system. Hence, for any reasonable model [Tf, Tg] = 0- O ne can calculate 
[Tf, T 2 9 ] employing the commutation relations for the Klein factors, the commutation relations Eq. i|16|) for the Bose 
fields and the Baker-Hausdorff formula. This results in [Tf, T|] = [T?JT$)+] = provided that Eq. (gSJ is satisfied. 
A different formulation of the same argument can be found in Ref. |7f Ref. discusses how Klein factors which 
ensure commutativity of tunneling operators can be derived using duality between weak quasiparticle tunneling and 
strong electron tunneling. 

The Klein factors serve as a manifestation of fractional statistics and are absent in the case of fermion tunneling, 
section III. A. The importance of Klein factors in quantum Hall systems with more than two edges has been empha- 
sized previousljiSiLiLiSiiS. In our problem, the Klein factors are necessary even though there are only two edges. 
Note that in the la) setup no Klein factors are needed. Indeed, the operators T\ and T 2 defined by Eq. (fT5f> com- 
mute. Commutativity of the quasiparticle tunneling operators is also ensured without Klein factors for the standard 
Aharonov-Bohm interferometer geometry, Fig. 3. 

In our calculations we will use the following representation of the Klein factors by 1/v x \ jv matrices: 



/ 1 ... \ 
1 ... 






V 1 







1 
... 1 

) 



/ v 
V 2 



K2 





V 1 







i) 1 !"- 2 

v 1 ^" 1 





(24) 



where ip = exp(— 2irvi). One can easily check that the above matrices satisfy the commutation relations l|23|) . The 
tunneling operators (|22|l with the Klein factors 12411 can be understood as products of a quasiparticle creation operator, 
an annihilation operator and a phase factor which includes the statistical phase accumulated during the tunneling 
event. 

The charge distribution on the edges is determined by the fields {4n(x)}, the total charge being determined by 
the zero modes of these fields. In the set-up depicted in Fig. la) the charge distribution completely describes all 
states in the Hilbert space on which the effective low-energy Hamiltonian acts. On the other hand, our discussion 
in section II shows that in the case of quasiparticle tunneling one needs to specify both the charge distribution and 
the effective statistical flux through the interferometer. Hence, the corresponding Hilbert space is the product of the 
space V^hargc on which the Bose operators <fii act and the space Vfiux on which the Klein factors act. As is clear from 



7 



the size of the matrices l|24|) the dimensionality of the latter space is dim Vfl ux = \ jv. This agrees with our discussion 
in section II where we found that the interferometer has \/v classes of states characterized by different probabilities 
of quasiparticle tunneling. 

Calculations based on the model H15|) with the tunneling operators (1221) confirm Eq. Q) of Section II. The Klein 
factors keep track of fractional statistics and are crucial for this result. In the absence of the Klein factors, i.e. for 
fractionally charged particles which do not obey fractional statistics, one gets qualitatively different results^. 

IV. ELECTRON TUNNELING 

We now study the electric current between the two edges. We consider the geometry depicted in Fig. la). The 
current operator 

7 1 1 

/= jfii = j.[H,Qi] = -^eM-ieVt/h)exp[i((t>i{0)/V^-MO)/V^)} 

+ l -^ L exp(ieVL/hv + 2t™$/$ ) cxp(-ieVt/h) exp[j(0i(L)/V^i - fa(L + a)/v^)] + h.c, (25) 

where Q\ is the total charge of the first edge and H the Hamiltonian. For our non-equilibrium problem we employ the 
Keldysh technique 2 ^. To this end we assume that the tunneling amplitudes Ti t 2 — at the moment of time t = — oo, 
and are subsequently turned on gradually. At t = — oo the system is in thermal equilibrium at temperature fcgT and 
chemical potential difference eV between the edges. The initial equilibrium state determines the bare Keldysh Green 
functions which will be used in the perturbative calculations below. 
The current at t = is 

I = Tr[pS{-oo,0)IS(0,-oo)], (26) 

where p is the initial density matrix and S(0, — oo) = Texp(— i J Hdt/h) the evolution operator. Expanding the latter 
to first order in the tunneling amplitudes one finds 

/= / ^(| ri |2 + |r 2 | 2 )exp(-^)[F(0,0,<)-F(0,0,-t)] 
n j_oo ft 

e f + °° ( ieVL ieVt 1 

+ / ddr 1 r*exp( 27ri$/$ ) exp( —)[F(-L,-(L + a),t) - F(L,L + a, -t)] + h.c. > , (27) 

ft J-oo { nv h J 

where 

F(b, c, t) = Tr[/5exp(i(/)i(x = b, t)/^/Tn) exp(-i0i(ir = 0, 0)/y/vi)] 

xTr[pexp(i0 2 (2; = c, t)/y/vZ) exp(- i 2 (x = 0,0)/^)]- (28) 

Let us first consider the zero temperature case. The correlation function is given by 21 

C ' l) = [S + t(t-b/v)}^ [8 + i(t-c/v)Yl»^ (29) 
where 8 is an infinitesimal positive constant, and r c is the ultra-violet cut-ofSS. With the above expression we find 

I = Io+h, (30) 

where 

/. = -— (iru' + lr.l-) ' ^ _ . (3D 
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2W /l/1+1/l/2 



exp(-27ri$/$ ) 



h 2 

1 d 1 /^- 1 exp(-^£) 1 d 1 '^- 1 exp(-H^) 

- 1)! (foiM-i '* =0 (0 + a/u) 1 /^ + (l/^ 2 - l)!dz 1 M-i |z =- a /" ?M 



1 x 



/ $ = " C 2 ( 1 )l/[2v 1 l+l/r2r/ 2 l : 

/l.C. 132) 



The current oscillates as function of the magnetic flux with period 3>o- At low voltages it follows the power law 

(see Appendix A; cf. Ref. fTo^) . Fig. 6 illustrates the J— curves for v\ = v-i = 1, = ^2 = 1/3 and v\ = 1, vp_ = 1/3. 
In the case when V\ = v%, similar expressions for other interferometer geometries were obtained in Refs. l5lll(l One 
can easily see that Eq. follows from Eqs. i|3QI32[l . 

Notice that only the difference a of the lengths of the edges enters the above expression for the current while 
the total length L of the first edge drops out. In the standard geometry of an Aharonov-Bohm interferometer— the 
flux-dependent contribution to the current depends on the total interferometer size and decreases with the system's 
size. Thus, quantum interference effects cannot be observed at large system's sizes. In the Mach-Zehnder geometry, 
quantum interference can be observed for L ^s> a since only the difference of the phases accumulated by the particles 
moving along two edges is important. 

The case of non-zero temperature is considered in Appendix B. The flux-dependent "interference" contribution J$ 
to the current vanishes at large aksT/ (hv). In the opposite limit of afc^T '/ '(hv) <C 1 as well as for ak B T ~ hv, the 
flux-independent contribution Iq and /$ are related by Eq. The linear conductance at low temperatures and low 
voltages eV <C k B T scales as 

G = I/V ~ [k B T] 1 / Vl+l / v - 2 - 2 . (34) 

V. QUASIPARTICLE TUNNELING 

We now consider the geometry depicted in Fig. lb). The current through the interferometer oscillates as function 
of the magnetic flux. In the first subsection below we determine the oscillation period. It turns out to be the same as 
in the case of electron tunneling and equals one flux quantum $o- Next, we use the perturbation theory to calculate 
the current as function of the voltage, temperature and the distance between the quantum point contacts. We confirm 
Eq. I|14|) . The dependence of the current l|14|) on the tunneling amplitudes T\ and T2 is nonanalytic. Such dependence 
cannot be obtained in any finite order of the perturbation theory. Thus, we have to sum up an infinite set of diagrams. 
In subsection V.B this is made for the simplest case of zero temperature and low voltage, eVa <^ hv, at the filling 
factor v = 1/3. The general case is considered in the final subsection. 

A. Period of Aharonov-Bohm oscillations 

The tunneling current operator can be found with the same method as in the previous section: 

/ = j t Qi = ~{H, Qi] = l -^El expt-ievVt/K) exp^(0i(O) - &(()))] 

+ ieK 2 vT 2 ^ levVL + 27r ^$ /$ Q ) expi-ievVt/K) exp[iv£(</>i (L) -fa(L + a))} + h.c. (35) 
n hv 

The average current is given by Eq. (j26(l . Since the tunneling amplitudes T; are small, we will employ perturbation 
theory. To lowest (second) nonzero order in Ti the resulting contributions are proportional to |Ti| 2 and | r2 1 2 • The 
cross-terms proportional to T*[T2 and TiTZ, vanish. Indeed, each term of the perturbative expansion is proportional 
to the product of the average of some function of the Bose-operators {(/>;} and the average of some product of Klein 
factors. The averages are determined by the initial density matrix p, Eq. (|26|l . The latter depends on the effective 
action Eq. (|15|l at t — —00. That action does not contain Klein factors. Hence, the density matrix p — p$p K , where 
p<f, and p K act on the spaces Charge and Vfl ux respectively, and p K is proportional to the unit matrix at any finite 
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temperature. Thus, the cross-terms are proportional to the expressions of the form Tr[/c;«;+], where I ^ m. Such 
traces are zero. This is readily seen from the following argument. For any two linear operators TiAB = TrBA. 
Hence, TrfftiK^] = Tt[«;^"ki]. At the same time it follows from Eq. I|2H[) that Tt[«;ikJ] = exp(27rz^i)Tr[K^"Ki]. Hence, 
Tr[KiK^] = 0. It follows that there are no cross-terms in the second order perturbation theory. 

The above result implies that the second-order terms are independent of the magnetic flux. The flux dependence of 
the current emerges only in higher orders of the perturbation theory. Each term of the perturbative expansion of Eq. 
(I26JI contains a product of Klein factors. Let the number of the Klein factors K? in a given term I a be n\ and the 
number of the factors Ki be nT ■ Each of the four numbers nf 1 indicates the power in which the respective coefficient 
Ti, T2, r* or enters in I a . Nonzero terms of the perturbation series describe the processes which do not change 
the edge charges. Hence, — n\ = — (n^ — n^)- 

We next show that vin^ — n^) 2 is an integer for any non- vanishing I a . Consider the trace W of the product of the 
Klein factors in the term I a . We can move all operators , k\ to the left of all operators rej, K2- This will produce 
a phase factor exp(i7). The trace can be represented as W = exp(i^)TiKiK 2 , where K\ is a product of ,K\ and 
K2 is a product of k\ , k-i . We know that 



TtK^ = TMTa^i- (36) 

At the same time, one can move all Klein factors k 2 5 ^2 in the product K1K2 to the left of the operator K\ using the 
commutation relations Eq. H23|) . This yields 



K-lK 2 = exp[27r^(n+ - n^) 2 ]K 2 Ki. (37) 

Eqs. I|36I37|) show that I a ~ TrKiK 2 7^ only if v{n J [ — n^) 2 is an integer. 

If v = l/[odd prime number], including experimentally relevant v — 1/3 and v = 1/5, the above result means that 

(rii —n>i) = — (nj — TI2) = q/v, where q is an integer. The term I a ~ T^ 2 (r|) n 2 ~ exp(27r^i$/$o["-^ — "2"])- Hence, 
it follows that the current is a periodic function of the magnetic flux with period <E>o- As shown in the following 
subsections the period is the same for any v — l/(2m +1). This can also be verified by a direct calculation of the 
trace ()36(l using (|24|l . Such periodicity agrees with the Byers-Yang theoremS&Si. If fractionally charged quasiparticles 
did not obey fractional statistics, i.e. if there were no Klein factors, the period would be $' = <&o/v (cf. Rcf. 16). 

Our effective hydrodynamic action I|15I22|) is applicable only for low temperatures and voltages. We are going 
to use the perturbation theory in I\2. It turns out that the lowest order contribution to the current scales as 
I r~j VT 2 2 max(eV, ksT) a , where a — 2v ~ 2 is negative. The perturbation theory can be used for the calculation of 
the tunneling current only when the tunneling current is much smaller than the current ve 2 V/h incoming from the 
sources. Thus, our calculations are valid provided that ri^max^V, fcsT)" -1 <C (jf^j ■ 

We will see that interference effects can be observed only for small enough a ~ hv /vaBxieV, ksT). This condition 
is similar to the restriction on the total interferometer size in Ref. 0- In our case the restriction on the total 
size L is weaker. As is clear from section II our analysis is based on the assumption that the time At between 
tunneling events exceeds the time L/v needed to a quasiparticle to travel from QPC1 to QPC2. Thus, L < vAt ~ 
vh 2 ( h 

2 t c ^max(eK/c£j T)t c 



B. i/ = l/3, T = 0, eVo<n« 

We will use the expansion of the current, Eq. (|26|l . in powers of the tunneling amplitudes Ti^. Only contributions 
proportional to even powers of the tunneling operators are non-zero. One might naively expect on the basis of power 
counting that the terms of order 2n in Ti a scale as V\Ti\P\T 2 \ 2n - p x [max(ey, T)] 2 ™^ -1 ) , if L - vh/max{eV, k B T). 
This is however not the case beyond the second perturbative order. In fact, the 4th order contribution (as well 
as higher orders) is infinite. To demonstrate this we use an analogy between the power expansion of the average 
current and the partition function of a Coulomb gas2&. Positively and negatively charged "particles" correspond to 
the tunneling events from edge lq to edge 2q and from edge 2q to edge lq respectively. The coordinates of the 
particles correspond to the times of the tunneling events. The charges can be located on both branches of the Keldysh 
contour, Fig. 7. The particles on the top and lower branches emerge from the expansions of 5(0, —00) and S(— 00, 0) 
in Eq. I|26() respectively. As discussed in Ref. 1261 attraction between oppositely charged particles binds them in pairs. 
The typical pair size is of order h/max(eV,kBT) at low temperatures and voltages^ (if L > hv/m&x(eV,kBT) then 
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pairs of size L/v are possible). When 



ri,2T c 

h 



x(eV,k B T)i 



< 1 and L < 



vh 2 ( 



2v-l 



the pairs 



are dilute and do not overlap. Hence, any term in the perturbation expansion reduces to the trace of a product of 
Klein factors times a product of the two-point correlation functions 



F[b = x% — x 2 ; c = [x% 



x 2 )(L + a)/L; t = tx — t 2 ] = Tr{/oexp[i\/W>i(x = Xx,tx)] exp[—iy/v(j)i(x 
Ti{pexp[iy/i'(j)2(x — xi{L + a)/L, tx)] exp[— i\fv<$> 2 {x = x 2 (L - 



= x 2 ,t 2 )]}x 
■a)/L,t 2 )]} (38) 



corresponding to the bound pairs. This product should be integrated over the times of all tunneling events, i.e. the 
positions of the charges. This integration can be separated into the product of the integrals over the dipole sizes 
(tx — t 2 ), Eq. (J3SJ, and the integrals over the dipole positions (tx + t 2 )/2. 

Fig. 7 illustrates the 4th order contribution to the current. One charge is located at t = and corresponds to the 
current operator in Eq. (|26[1 . It forms a dipole with an opposite charge at the point r ~ h/m.ax(eV, kgT). Two 
more mutually opposite charges are located at the points tx ±n/2, where n is the dipole size. The charges can reside 
on the same or different branches of the Keldysh contour as shown in Figs. 7a), 7b) and 7c). The integral over tx 
diverges. In the absence of the Klein factors the infinite integrals corresponding to the configurations Fig. 7a), 7b) 
and 7c) cancel but this is no longer the case when the Klein factors are included. Hence, the 4th order contribution 
is infinite. The same argument applies to higher-order contributions. Certainly, the sum of all perturbative orders 
must be finite but as is clear from section II it is not an analytic function of I\2- Below we develop a method to sum 
up all orders of the perturbation theory. 

We begin with the simplest case when v = 1/3; T = 0; eVa, eVL <C %v. The same approach will be applied to the 
general situation in section V.C. The first condition simplifies the structure of the Klein factors l|24|l which become 
3x3 matrices at filling factor 1/3. The second condition allows us to use a simpler zero-temperature expression for 
the correlation function 



F(b,c,ti — t 2 ) = Tr[pexp(i v / ^^i(x = b, t\)) exp(— i^/v4>i{x — 0,t 2 ))]x 



Tr[pexp(z^^ 2 (a; = c, tx)) cxp(-iy^(j) 2 (x = 0, t 2 ))] = 



[5 + i(ti - *a - b/v)) v [S + i{ti - *a - c/v)Y 



(39) 



The third condition makes it possible to neglect the distances L and (L + a) between the point contacts. Indeed, the 
tunneling operator Tg =: T 2 k 2 exp(iy^[(/>i(L, t) — <fi 2 (L + a, t)]) : +h.c, Eq. (|22|l . can be rewritten as 



d k (/)i 



T| =: r 2K2 exp ( V^l(0,*) - <M0, t)\ + 



fe=i 



k dkh (L + a) k }/k\] :+h.c. 



: r 2 K 2 exp(i v ^[</'i(0,i) - «^(0,t)]) {l + iV^ 



dx k 
d<j) 2 



U<pi U>p 2 



dx 



dx 



+h.c. 



(40) 



where the ellipses denote higher order gradients of </>/ and the colons denote normal ordering. After the substitution 
of Eq. (|40|) in the perturbative expansion of the current one can compare the contributions from the terms containing 
derivatives of <pi with the contributions from the terms which do not contain such derivatives. Power counting 
shows that the terms with derivatives are suppressed by the factors of order (LeV/hv) k ; ([L + a]eV/hv) k and hence 
can be neglected. This conclusion agrees with the results of section V.C for arbitrary voltages, temperatures and 
interferometer sizes. 

Thus, we can use an effective single impurity model. The tunneling operator O = {Tf +T 2 9 ), Eqs. (|22I24II . assumes 
the form 



= 0-+0+; 0_=0+=/t_A_, 



(41) 



where 



A- = A\ = exp(V^i(0,i) - 0a(O,t)]) 



(42) 



and the (non-unitary) operator k_ is 
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/ Ci \ 

k_ = k ] + = Tiki +r 2 exp(27ri/i$/$o)K2 = C 2 (43) 

V^3 / 

with Ci = Ti + r 2 exp(27rw$/$o)'i/', C 2 = I\ + r 2 exp(2irvi$/$ )ip 2 , C 3 = Ti + T 2 exp(27r^i$/$ ) and ip = 
exp(— 27ri/3). In what follows we will denote the basis vectors (1, 0, 0), (0, 1, 0) and (0, 0, 1) as (1|, (2| and (3| 
respectively. Since we neglect L and [L + a) in the rest of this section, we can use the correlation function l|39|l in the 
simplest limit b, c — 0. 

The contribution I 2 m of order 2N to the current corresponds to the charge configuration with N dipoles. One dipole 
of size ro is located at t = 0. The remaining dipoles of sizes t\, . . . , tn-i are located at > t\ > t 2 > . . . > tisf—i. 
One finds 

hN = - ^2 ^2 ^2 cr[0] 2N I dtx I dt 2 ... dtN-i / dr / dn . . . dT A r_ 1 & II/ c &+6~ 

6 =±i 6±=±i <r[0]=± J -°° J -°° J ~°° J -°° 

Tr [pT c O a[0] (t = 0)O_ CT[0] (r ; 6 )O_(ii + -; -y,b+)... 0_(i w _i + ^yi; ft^O+^-i - 6+„ 1 )](,44) 

where T c denotes time-ordering along the Keldysh contour, bo, bt, b k = +1 correspond to the bottom branch of the 
Keldysh contour and bo,b~£ ,b k = —1 correspond to the top branch. 
The trace in Eq. 144|) can be factorized as 



Tr[...] = n A n«, (45) 



where ±1,4 stays for the trace of a product of operators A± , and II K denotes the trace of a time-ordered product of the 
Klein factors. The former trace can be further factorized into a product of the two-point correlation functions 1391) 
corresponding to each dipole. We next simplify the expression for iI K , 



n K = Tr 



PkTcK^q] {t = 0)k_ <7 [o](t ;6o)k-(*i + y^i - • • •«-(%-! + -y-^-; b^.^K+^iV-i ~ -y^jv-i) 



(46) 

Using the matrix elements of the Klein factors the above equation can be rewritten as 



TT \ " \ " Tl t b\r>l t b t b\ r>( t b t b \r>/ t b t b \ 

11k = 3 J(ei,e 1 )i?(e 1 ,e 1 ;e 2 ,e 2 ) . . . B{e N _ 2 ,e N _ 2 ; e N _ 1 ,e N _ 1 )B{e N _ 1 ,e N _ 1 ;e N = e N ,e N = e N ), 

(47) 



3 

ejv-1,2,3 e ?/ 1, =1.2,3;k=l,...JV-l 



where the factor 1/3 comes from the condition Tvp K = 1; I(e\, e\) = (ei\T c n a [ Q }(t = 0)k_ <t [o](to; &o)|e*) an d 



B(e<,4;4 +1 ,4 +1 ) = (e^ +1 |e 6 fe )(e*|T cK _(t fc + r fc /2;-l) K+ (t fc -r fe /2;-l)|e* +1 ) for 6+ = 6 fe T = -1 

B(e*,e|;e^ +1 ,e| +1 ) = (e£ +1 |T cK _(t fe + r fc /2; +l)/s+(t fc - r fe /2; +l)|e^)(e t fc | e *. +1 ) for 6+ = 6" = +1 

B(e*,e|;ei +1 , e b fc+1 ) = - r fe /2; +l)|e 6 fc )<4 |/c_(t fc + t*/2; -l)|e£ +1 ) for 6+ = +l;6^ = -l 

B(4,e*;4+i>4+i) = (4+ilM** + r fc /2; +l)|e b fe ) (4M*fc - r fc /2; -l)|e* +1 ) for 6+ = -1;6^ = +1 (48) 

It is clear from Eqs. I|48l) that for equal e k+1 — e b k+1 — e k+ i, the expression B(e k , e b k ; ek+i, ek+i) is nonzero only if 
e\ = e\. Hence, Eq. Q47JI can be represented as a matrix product 



II K = ^2 ^(ei)-B(ei, e 2 ) . . . B(e N - 2 , e N -i)B{e N -i,e N )p{e N ) 1 (49) 

efc 

where p(e.N) — 1/3, -B(efc, efc+i) = B(e k ,ek',ek+i,ek+i) and /(ei) = I(ei,ei). After the substitution of Eq. I|49|) 
in Eqs. I|45l) and l|44|). each matrix element .B(e/., e&_|_i) multiplies by the correlation function i 2 (T c b k A_(t k + 
Tfe/2; &^T)6 ; tA + (ifc — Tfc/2; &£)). The product should be integrated over dru- One finally obtains 
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I 2 n = / dti dt 2 

J — oo J — oo 



(50) 



where 



(p| = (1/3, 1/3, 1/3); 



(51) 



h ^ 

fe=l,2,3 



(k\K+K-\k) 



exp(ipeVt/ti)Tc 
(S + it) 2v 



(k\K-K+\k) 



exp(iveVt/h)i 



{6 - it) 2 " 



(k\ 



^(-)- M|CsMCl|2 , m 



(52) 



D 



V3r(l/3)r c 2/3 feV 



-1/3 



|c 3 | 2 





-|Ci| 2 
\Ci\ 2 



|C 2 I 2 



(53) 



The total current 



I = ^2 L 2N = (I\cxp(iD)\p), (54) 

where i — +oo is the length of the Keldysh contour. 

All elements of the matrix D are real, all diagonal elements are negative, all nondiagonal elements are positive or 
zero and the sum of the elements in each column is zero. The Rohrbach theorem^ applies to such matrices. Since 
(p\D = 0, is an eigenvalue of the matrix D. According to the Rohrbach theorem, this eigenvalue is non-degenerate 
and the real parts of all other eigenvalues are negative. This allows for a simple calculation of exp(£D). Let S be such 
a matrix that D = SDS^ 1 assumes the Jordan normal form. Without the loss of generality we can assume that the 

first column and the first string of the matrix D are zero. Then the first string of the matrix S can be chosen in the 
form (1, 1, 1). Thus, the matrix exponent 












1 




• 












• 









v o 








exp(LD) = S~ l j n (! i) I S = .S-' j li n (i j ir,r,j 

Hence, 



/=(/|5- 1 |l) = (/|5r 1 ), (56) 

where IS'-f 1 ) denotes the first column of the matrix S^ 1 . 

In order to complete our calculation we have to determine the components f\, f 2 , and of the vector IS*^ 1 ). From 
the condition SS^ 1 = E one finds that 

h + h + fa = 1. (57) 
We also know that IS^ 1 ) is the eigenvector of D with zero eigenvalue, i.e. £)|S'^ 1 ) = 0. Hence, for each k, 



\C k -i?h = \C k \ 2 f k+u (58) 

where we use the convention 3 + 1 = 1. The solution of the above equation is ff. — a/\Ck-i\ 2 , where 
a = i/ \ c 1 \i+i/ \ c 2 \' i +i/ \ C3\' i can ^ e found from Eq. JS7J. Finally, 
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-er(l/3)r c 2/3 feVY 



1/3 o 

i T i T i ■ (59) 



This expression is equal to the harmonic average of the three tunneling currents in three systems with a single quantum 
point contact with the tunneling amplitude C\, C 2 and C3 respectively. In other words, Eq. 1(59(1 is equivalent to P|. 
Using the result of Appendix C, Eq. (|5$f can be represented as 

' = ~ er< vit?'" (S") " 3 [|ri16 + |r2l ° + 2|rir2 ' 3 cos(3ao + 2,r * / * o)1 |fl} 



.v (60) 

where a = arg[r 2 /Ti]. This allows one to easily verify Eq. (0) in the case u = 1/3, T = 0, eUL, eUa <C 

C. Tunneling current for arbitrary filling factors, temperatures and voltages 

The calculations follow the same route as in the previous subsection. The contribution to the current of order 27V, 
I2N, expresses via the trace II of the time-ordered product of 2N tunneling operators which form N dipoles. The trace 
must be integrated over the size of each dipole and the positions of N — 1 dipoles, the remaining dipole being located 
at t = 0. The trace II factorizes as the product II = n K Il0, where 1I K stays for the trace of the product of the Klein 
factors m, K2, k+ and kJ; II^ stays for the trace of the product of the operators exp(±iy/v[(j)i(x — 0) — 2 (x = 0)]) 
and exp(±i v ^[^ 1 (L) — <fi 2 (L + a)]). The latter trace factorizes in the product of the two-point correlation functions 21 
corresponding to each dipole: 

F(b,c,ti - t 2 ) = Ti\pex^{iy/v<j>i{x = b, ti)) exp(—i^^i(x = 0, i 2 ))]Tr[/5exp(zV^'/ , 2(a; = c, ti)) exp(—i^(j>2(x = 0,i 2 ))] 



nkBTT c /h 



sin(irk B T[6 + i{h - t 2 - b/v)]/h) 



nk B TT c /h 



sin{irk B T[S + i(ti -t 2 - c/v)]/h) 



The former trace II K can be represented in the form similar to Eq. 149(1 . where B(e k , e k +i) express via matrix elements 
of Ki and K2- Next, one can rewrite I 2 n in the form 1(50(1 with modified definitions of \p), (I\ and D. \p) and \I) are 
now 1/^-dimensional vectors; D is a matrix of size \jv x 1/za Similar to the previous subsection, in order to obtain the 
matrix elements of D, one has to multiply B(eu, ek+i) by the correlation function 1(61(1 describing the dipole located 
at t = tk and then integrate the product over the dipole size. One finds 

(p\ = {v, v,..., v); (62) 



(I\e k ) = ev{I k{V) - Jfe(-v)); (63) 
(ei\D\e k ) = 5 k ,i(Ik(-v) + h(v)) - h,i-ih{-v) - h,i+ih(v), (64) 

where 

h{-v) = i [d r i| 2 + |r 2 | 2 )j(-^; 0) + T.THrfji-V; -a) + j(-V; a)} ; (65) 

n 



h ( v) = [d r il 2 + |r 2 | 2 )j(^; 0) + r^r^-'jiv; a) + rir^jiv-, -a)] , (66) 

where ip = exp(— 2-niv), the phase factor exp(27rii/"I>/$o + ievVL/fiv) should be included in r 2 and 
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r+oc 

j(U,0) = - I dtexp(iueUt/K)F{0,0,t); 



/-t-oo 
dtexp(iveUt/h)F(±L, ±{L + a),t) (67) 
-OO 



The total current is given by Eq. (|f>4|> of the previous section with the above definitions of D and \p). 
One can easily check that j(U,x) — j*(U, —x). Hence, Iy±v) are reai - Appendix D shows that 



/ veV\ 

h(-v) = exp I ~^7p J I (k+i)(v) = ll(k+i)(v) (68) 

that is equivalent to the detailed balance condition JSJ. A comparison with the geometry Fig. 3, Ref. 0, shows that 
the same integral which expresses the current ^Fig.3 m the geometry Fig. 3 also equals the expression 



-/ (fe -i)(-v) = (i-7)4(i/) =4ig. 3 /K> = — ^-|r efr | — ^ — J Y{2v) 2fc^T' (69) 

where 



|r ofr | 2 = 2tt-^ V r m I7 exp(i[27ri/$/$ - euVa/(2hv) - 2ir(k - l)H[<5 m , 2 - tf J)2 ])x 

exp(— virak B T /[Tiv]) ( eyLp(ievVa/[2Tiv\)F (u, v - ievV/2itk B T; 1 - ieuV/2irk B T; e -^k B T a /[hv]} 
smh(eVv/2k B T) m \ I> + iez^/27rfc B T)r(l - ievV/2irk B T) 

F being the hypergeometric function^. 

Since the particle current flows from the edge with the higher chemical potential to the edge with the lower potential, 
we know the sign of Jptg 3 = dQx/dt > (electron charge e < 0). Hence, Ik(±v) < 0- One can also see that the 
sum of the elements in any column of the matrix D, Eq. (|64fl . is zero. Thus, the Rohrbach theorem^ applies again 
and hence exp(il)) = S~ 1 MS, where the matrix S reduces D to the Jordan normal form, the matrix M has only 
one non-zero element in its upper left corner and this element is equal to 1 just like in the preceding subsection. All 
elements of the first string of the matrix S are equal to 1. This allows us to obtain the expression for the current in 
the form (|5fc>|) following exactly the same steps as in section V.B. 

The only thing left is the calculation of the components f r of the vector l^j -1 ), i.e. the first column of the matrix 
From the condition SS^ 1 = E one finds 



r=l 

From the condition DIS^ 1 } = one finds 



h(v)fk - I(k-i)(-v)fk-i — I(k+i)(v)fk+i - h(-v)fk- (72) 
At the same time the current (|63I5GH is 

I = vey^Hv) ~ h{-v))fk (73) 
k 

The system of equations (|71I73|) together with the detailed balance condition l|t)8|) is equivalent to the system H8I1U|) 
of section II with Ik(±v) playing the role of the transition probabilities p. The components f r of IS 1 ^ 1 ) have the 
physical meaning of the distribution function. Thus, we can directly use the solution l|14|) : 
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1=^^- , (74) 

where the current I' k = Ipf^^, Eq. I|()9|) . equals the tunneling current in the geometry Fig. 3 with the tunneling 
amplitudes Ti and T 2 exp(z[27w<E>/<E> - evVa/{2Tiv) - 2n(k - l)u]). 

In the limit T — 0; aeV <C hv the above result reduces to a simple generalization of Eq. ijfJUJ): 



-2ei/siii(27ri/)r(l - 2v)t* v ( eVv\ |r x | 2 - |r 2 | 2 



|ri| 2 /"-|r 2 |2/" 



|ri| 2/ " + |r 2 | 2 /" + 2|r 1 | 1 /^|r 2 | 1 /^ cos(27r$/$ + a /v) 

(75) 



If |Ti| ^ |r 2 | the current never vanishes. At |Ti | = |T 2 | a "resonance" is reached when $ = [n + 1/2 — ao/(27r^)]$ 
and J = 0. 



VI. CONCLUSIONS 



We have found that the tunneling current through the electronic Mach-Zehnder interferometer is a periodic function 
of the magnetic flux with period $o. This result is valid both in the weak quasiparticle tunneling regime and in 
the weak electron tunneling regime. The relations between the flux-dependent and flux-independent components 
of the current, 7$ and Iq, are different in these respective regimes. In the electron tunneling case, I$(ri,r 2 ) ~ 
[Jo(ri,r 2 ) — Jo(ri,0)] 1 / 2 at low voltages and temperatures, T 2 -C IV In the quasiparticle tunneling case the flux- 
dependent contribution scales as i$(IVr 2 ) ~ [io(IVr 2 ) — Io(Ti,0)] b , where b > 1. The exponent in this power law 
contains information about quasiparticle statistics since the exponent derives from the algebra of the Klein factors. 
Thus, the Mach-Zehnder interferometer can be used to probe fractional statistics. Recently an interference pattern has 
been observed experimentally for the integer quantum Hall case in a Mach-Zehnder interferometer—. Higher magnetic 
fields would allow an investigation of a fractional quantum Hall liquid. 
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APPENDIX A: ELECTRON TUNNELING AT LOW VOLTAGES AND ZERO TEMPERATURE 



In this appendix we verify that the current, Eqs. (|30I32() . satisfies the asymptotics <|33ll at V — > 0. The flux- 
independent contribution to the current 1|31|) scales as V 91+92-1 : where g\ = \jv\ and g 2 = l/v%. Hence, we need to 
check that the flux-dependent contribution J$ does not exceed constT^ 91+ff2_1 at low voltages. Expanding exp(— -^p) 
in lj%2l in powers of ieVz/Ti one finds that the coefficient before V k in the Taylor expansion of 7$ equals 



Sfe = const 



d 92 - 



(, 9 i - 1)! dz 91 ' 1 1 (z + a/v) 9 * (g 2 - 1)! dz 



yk-gi 



(Al) 



We want to show that Sk = for all k < g± + g% — 1 . 
Let us consider separately k > g\ — 1 and k < gi — 1. 

1) k > gi — 1. This inequality can be satisfied for g 2 > 1 only since k < g\ + g% — 1 
\ z =-2.z k - gi = since g 2 - 1 > k - g x > 0. Since £^\ z = a z k = at p < k, one finds £^ 1 



9l-i-dzP~ dzm- 1 -P (z+a/v) g 2 

.91 - 1. 

2) At k < gi — 1 one finds 



In this case 
- z " - 

r|z=0 ( z+a /v)B2 — 

= 0, where C b a denote binomial coefficients. Thus, Sf. = for gi + g 2 — 1 > k > 
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Sk = const 



i rk fcK-i) 91 " 1 "" (gi + g 2 -fc-2)! 

(<7i-l)! 9l ~ 1 (a/v)9i+92-k-i (. g2 -l)! 



1 



(.92 - 1)! 



(-a/v) 



k- Ql - a2 +i (9i + 92 - k - 2)\ 



(91 



k-l)\ 



Thus, sfe = for all k < g\ + g% — 1. 



= 0. 
(A2) 



APPENDIX B: PERTURBATION THEORY AT FINITE TEMPERATURE 

We want to calculate the current, Eq. lj2*7)l. at a finite temperature, k B T > 0. The correlation function, Eq. H28JI. 

-.21 



F(b,c,t) = 



Trk B TT c /h 



sm(Trk B T[5 + i(t - b/v)]/h) 



TtksTTc/H 



sm(nT[6 + i(t - c/v)]/h) 



(Bl) 



Since = gi and 1/^2 = 92 are odd integers, F(—L, — (L + a),t) = F(L, L + a, — t) on the real axis except in the 
vicinity of the two real poles t = L/v and t = (L + a)/v. F(0, 0, t) = F(0, 0, —t) except near t = 0. Hence, the integral, 
Eq. I|27l) . reduces to the sum of the contour integrals along small circles around the real poles. A straightforward 
calculation yields 



/ = 



"" /7rfcsTTcV "(|r 1 | 2 + |r 2 | 2 ) Jo + 



er c / Trk B TT c \ 9 1 



[r x r^ exp(-27ri$/$ ) J$ + c. 



(B2) 



where g = gi + g2 = l/^i + 1/^2, and 



Jo 



2tt 



i 9+1 (g~ 1)1 dzs- 1 



\z=0- 



exp(ieV z / irksT) z 9 
sinh 9 z 



(B3) 



J$ = 



2tt 



2=0 



?9i 



2tt 



d 92 



(sinh 2;) 91 (sinh[z — 



ak B T - 

hv 



)92 i9+ 1 (g 2 -l)\dzS 



hv ) 



(sinh z) gi (sinh[z 



2fe„T- 



92 



(B4) 

Note that the temperature enters the above expression in the combination ak B T but not in the combination Lk B T. 
This can be understood from the picture of non-interacting electrons. For non-interacting particles the current is the 
sum of independent contributions from different electron energies. Each contribution depends on the phase difference 
between the two paths connecting the point contacts but not the total phase accumulated on any of those paths. The 
former is proportional to a, the latter to L. 

One can easily extract the asymptotical behavior of the linear conductance G — ^ at low temperatures (|34H from 
equations (|B3IB4I) . In the Loran expansion of G in powers of a, each term of order a™ is proportional to j ,1 / v i+ 1 / 1/ 2-2+n 
since a enters the expression for the current in the combination (ak B T/hv) only. In the limit a — > the conductance 
must remain finite. Indeed, at L = this limit corresponds to a problem with a single tunneling contact. Hence, only 
positive and zero powers of a are present in the Loran expansion and the leading contribution to the temperature 
dependence of the conductance scales as T 1 ^ Ul+1 ^ 2 ~ 2 . 



APPENDIX C: EFFECTIVE SINGLE IMPURITY MODEL 



In the effective single impurity model, section V.B, the current is given by Eq. (|59H . The purpose of this appendix 
consists in the calculation of the coefficient in Eq. H59JI . 
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where C k = I\ + T 2 exp(27ri$/[iV$ ] - 2mk/N), Eq. |@3J|, iV = 1/za Without the loss of generality we can assume 
that Ti = \Ti \ is real. Let T 2 = |r 2 | exp(zao). In the rest of the appendix we will use the notation ji — |r 2 |, 7 2 = |T 2 | 
and <j> = 27r$/ [7V$ ] + ao- U can be represented as the ratio of two polynomials of 71 and 7 2 : 

TJ = N HfelTi + 72 exp(i^) + 2mk/N)\ 2 

E fe n; #fc |7i+72exp^ + 2^/iV)P' 1 J 

where the prime after the product sign means that the term with I = k is not included in the product. The nominator 
in Eq. I|C2|) equals |P( 7 2 exp(icf>))\ 2 , where the polynomial 

P(z) = ILJ71 + zexp(27rifc/7V)] (C3) 

All roots of the polynomial i|C3(l coincide with the roots of the polynomial 7^ + z N . Hence, from the basic theorem 
of algebra 

P(z) = 7 f + z N (C4) 

and the nominator is 

|P( 72 exp(^))| 2 = (li + 12 exp(iiV»)( 7 f + 7 * exp(-ti\ty)) = 7l 2W + ^ N + 2 7 f 7 * cosNcf,. (C5) 
The denominator in Eq. I|C2|I can be represented as 

d = E gfr^gg) 2 = 2 I £ (_l), 7 f -i-p^ exp( lp + 2^/iV)| 2 = 

7i + 7 2 exp(z<p + zmk/N) ^-^ ^— ' 



k p=0 



k 

N-l JV-1 iV-1 

^ ^ li~ 1 ' P ll{~l) p e^{ip4) + 2nkpi/N) x 7f- 1 - r 7^(-l) r exp(- 2 r ( /.-2^rz/iV). (C6) 

fe=0 p=0 r=0 

The sum J2k=o exp(27r/c(p — r)i/N) = 0, if p — r 7^ niV. With the help of this property, Eq. (|C6|1 reduces to 

N ~ 1 „,2iV -,2JV 

rf = ^E T?^- 1 '- 2 ^ = ^^V^- (C7) 
Finally, the combination of Eqs. ljC5f> and (|C7|) yields 

2 2 

rr 7l ~ 72 r 2N . 27V . r, JV JV M il //-ion 

17 = 2jV_ 2N Ul +72 + 2 7l 72 cos N 4>\ (C8) 

7i 72 

APPENDIX D: DETAILED BALANCE 

Here we derive Eq. (|68|l . We need to show that 

j(-V; 0) = 7 j(V; 0);j(-V; ±0) = 7 j(F; T«) (Dl) 

This is equivalent to the equation 

r+°° ievVt f +oc ievVt 

J dtF(-b,-c,t)exp( — )= 7 y dtF(&,c,i)exp(_-— ), ( D 2) 
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where F is given by Eq. i|61|) . Both integrals are taken over contour a) which goes below the real axis in Fig. 8. We 
can change the sign of t in the first integral in Eq. I|D2|) . This also changes the integration contour into contour b), 
Fig. 8 (infinitesimal S is positive). The integral over contour b) is equal to the integral over contour c). The latter 
integral equals 7 J_°° dtF(b, c, t) exp( iei ^ vt ) indeed. 
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FIG. 1: Schematic picture of the Mach-Zehnder interferometer. S and D denote sources and drains. Arrows indicate the 
propagation direction of the chiral edge modes. The source and drain voltages are Vsi = V, Vm = Vd2 = Vs2 = 0. a) 
Electrons tunnel between the puddle defined by edge 2e and the QHE strip (edge le) b) Fractionally charged quasiparticles 
tunnel between edges lq and 2q. In both cases the difference of the chemical potentials equals the voltage drop between SI 
and S2. 




FIG. 2: Quasiparticle propagation and tunneling in the Mach-Zehnder interferometer. Once quasiparticle q\ a with its attached 
flux tube tunnels through QPC2 it arrives in D2 (qit) and the flux tube is released within the MZ hole. The flux tube 
affects the statistical phase accumulated by the next tunneling quasiparticle. A quasiparticle emitted from S2 will have its flux 
tube folded (525) and will not affect the statistical phase accumulated by other quasiparticles. 
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FIG. 3: The Aharonov-Bohm interferometer from Ref. 0- 
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FIG. 4: The Aharonov-Bohm interferometer considered in Ref. 7 
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FIG. 5: Mach-Zehnder interferometer can be modeled as two chiral Luttinger liquids with the same propagation direction. 
Dashed lines represent quantum point contacts. In the case corresponding to Fig. lb) quasiparticle tunneling is allowed while 
in the case that corresponds to Fig. la) only electron tunneling is allowed. 
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FIG. 6: Current I(V,$ = 0) (thin line) and visibility R = [max*I(V,<&) - min* /(V, $)]/ [max*I(V, $) +min*7(V,$)] = 
max/cj//o (thick line) as functions of V/Vo = eVa/(%v). T he cu rrent is normalized to its value at V = 25Vo; Ti = T2; a) 
H = V2 = 1; b) za = 1, !/ 2 = 1/3; c) 1/j = v 2 = 1/3. [cf. Refs. iBllfi ] 
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FIG. 7: The 4th order contributions in the Keldysh technique can be represented in terms of two dipoles one of which is located 
at t = 0. The second dipole can include a) two charges from the top branch of the contour; b) two charges from the bottom 
branch; or c) one charge from the top branch and one from the bottom branch. 
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FIG. 8: Three integration contours in the complex plane which we use in the derivation of the detailed balance condition. 



